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Abstract: 

We have analysed the loffe-time distribution of quarks in virtual photons 
using Operator Product Expansion of the correlation function that deter- 
mines the matrix element of the corresponding quark string operator. The 
distribution for a transversally polarised photon admits a spectral represen- 
tation which can be continued to the on-shell region = 0. The resulting 
model loffe-time distribution turns out to be larger than parametrisations 
of the available F2 data. This result is linked to the slope of the quark 
distribution at the origin, which is also too large. 
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Many years ago Witten [|1| showed that the Operator Product Expansion 
(OPE) can be apphed to the description of photon-photon scattering with 
the target photon on or near mass shell while the probing photon has a 
large Euclidean momentum. According to Witten's analysis the structure 
function of a transversal polarised photon contains a point-like part which 
grows logarithmically as a function of . It arises due to the contribution 
of highly virtual quarks in the quark loop and can be fully predicted in 
perturbative QCD. The contribution of low-virtuality quarks describes the 
interaction of the probe with the hadronic component of the photon, which 
is subdominant at asymptptically large . 

It is necessary to understand this hardronic contribution to the photon 
structure function better because all experiments probe the photon struc- 
ture at values, for which its hadronic part is non-neglible. Because it 
has not yet been possible to calculate it from first principles. Vector Me- 
son Dominance (VDM) type models are used to estimate its magnitude. To 
link photon properties more directly to QCD, Gorsky et al. calculated 
the hadronic part of the photon structure function applying OPE to the 
photon-photon forward Compton amplitude, which resulted in an expansion 
in inverse powers of the target photon virtuality. This expansion admits a 
spectral representation which allowed an analytic continuation to the quasi 
real target photon regime ~ 0. 

However, the problem with parton distributions in momnetum space is 
that for a given value of the momentum fraction u they get contributions from 
both small and large longitudinal distances. To circumvent this problem the 
authors of 0, Q proposed an alternative representation, which was already 
applied to the quark distribution of the nucleon 0. It turned out that in 
complete duality to twist-2 parton distributions in momentu space, twist-2 
parton distributions in coordinate space can be introduced. The coordinate 
space variable, called loffe-time, measures the longitudinal light-cone distance 
between the points where the hard probe was absorbed respectively emitted 
by the target. Derivatives of loffe-time distributions at the origin are given 
by moments of some coresponding structure functions, or equivalently § 
by matrix elements of twist-2 operators of increasing dimension. From the 
theoretical point of view, this connection makes loffe-time distributions much 
easier to analyse than parton distributions in momentum space. 

In the present paper we concentrate on quark distributions q^{u) and 
q'^iu) in a photon with momentum p and polarization A. Following reference 
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[^] we introduce the loffe-time distribution of quarks inside the photon Q^{z), 
which is given by the photon matrix element of the quark string operator 
6(A) = ^(A)^ [A; -A]^(-A), 

2i{p-n)Q\z) = {pX\d{A)\pX) 

r-l 



2{p ■ n) 



du 



q\u) 



exp yiuz) 



q'^{u) exp 



-lUZ] 



where for notational convenience we have introduced two parallel light-like 
vectors A^ and satisfying A^ = = n • A = 0. Our convention is such 
that for any vector a, n ■ a = a+. In equation z = 2p ■ A and [A; —A] 
denotes the path-ordered exponential which ensures gauge invariance of the 
quark distribution. As we have to respect both QED and QCD p gauge 
symmetry [A; —A] becomes the product of colour and electromagnetic parts 



[A; - A] = Fexp 



d^A ■ A{^A) 



Fexp 



— ze 



dT]A ■ A{t]A) 



(2) 



where v4^(x) 



is the gluon field and v4^(x) denotes the photon field. In the 



following, however, we will adopt the Schwinger gauge x^A^ 



X 



for the 



electromagnetic field, such that the second factor becomes one, and never 
explicitely enters the calculations. 

A Taylor expansion of 0(A) in A^ results in local operators of twist 2. 
The first term in this expansion is given by the traceless part of the quark 
energy-momentum tensor and the matrix element of its QCD part describes 
the longitudinal momentum fraction carried by quarks in the target. The 
normalization of loffe-time distributions is such that this matrix element is 
given just by the derivative of Q^{z) at z = 0. As the photon is even under 
charge conjugation (|l]) becomes 



{pX I 6(A) I pX) = 2i{p' 



n] 



du 



q\u) 



u 



simuz] 



(3) 



which clearly shows that the loffe-time z is the Fourier-conjugate to the usual 
Bjorken variable u. 

Using the LSZ reduction formula we calculate the photon matrix element 
(|1|) as a three-point correlation between two electromagnetic currents and 
the operator 6(?/; A) centered at a point y: 



d{y- A) = ^!{y + A)i[y + A-y- A]m{y - A) 



(4) 
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Denoting by M^y{p^ k) the three-point correlation function 

M^,,{p,k) = i'e' J d^xd%exp{ipx + iky){0\T{j^{x)d{y;A)j,{0))\0), (5) 
we obtain 

{pX I Oiy; A) | pX) = e;(p, X)M,,{p, k = 0)e,{p, A) (6) 

We remind the reader that (^) holds only in the Schwinger gauge for the 
photon field; in general some extra terms arise. 

With the Ward identity derived in the appendix of we can write 
M^y{p, k) as a sum involving two- and three point correlators: 

M^,(p, A;) = 

j d^x exp {ip ■ x) j d^y exp {ik ■ y) j^— J (0 | T^jf^{x)'^{y + A) 
x[y + A;y + ^A]gYA'^Gp^{y + ^A)[y + CA;y- A]^il{y - A)j.(O)} | 0) 
d^x exp {ip • x)- exp {ik • a:)(0 | T{j^{x] x — 2A)j^(0)) | 0) 

K, ' X 

- (0 I T{j,{x)3,{2A- 0)) I 0) + (A ^ -A)] , 

(7) 

where j^{x,y) is the point-splitted electromagnetic current operator, 

j^{x,y) = ^{x)^^[x;y]<i/{y). (8) 

Note that the momentum k^ appears in the denominators in and therefore 
it can be taken to zero only at the end of calculation after all potentially 
divergent terms canceled. Expression is formally equivalent to (||) but 
allows to apply directly the techniques developed previoulsy in 

In the following we distinguish between loffe-time distributions in transver- 
sally (T) and longitudinally (L) polarised photons, defined by 

^ 5] (pA I 6(A) I pX) = 2i{p ■ n)Q^(z) 



2 

^ A=l,2 



{pL I 0(A) I pL) = 2i{p ■ n)g^(z) (9) 
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We introduce the vector 



P.=P.-^^n,, (10) 



such that = 0. The polarisation tensor can than be written in terms of 
transverse and longitudinal parts, d^i, = — d^^: 

dlu = H <b> A)e^(p, A) = -g^^ + ^^^^ ^ ^"^'^ 

A=l,2 P ' 

< = -'>^^>^^P^^^ = -P'^2- (11) 

Our main task is using OPE to calculate the three-point correlation func- 
tion for the case of massless quarks up to the leading non-perturbative 
corrections of dimension 4, 

Q^'\z) = Ql'\^) + j^,Ql\z) (12) 

where {—p^) is the virtuality of the target photon and A its polarisation. 

We shall first consider the situation for transverse polarisation. The per- 
turbative contribution Ql{z) to the matrix element (|^) arises from the last 
two terms in (0). The explicit calculation gives: 



QI{z) = — [ du sin(M-2;) 

71" „ ^ J 



m2 

(m^ + u^) log (At) - {u^ + u^) log (uu) - 1 

(13) 

where is the ultraviolet cut-off corresponding to the normalisation point 
of the operator O, and a = ^ is the elctromagnetic coupling constant. 
Assuming a relatively low normalisation point ~ a few GeV^, the sum 
over active quark flavours runs over u,d and s quarks and = |. Note 
that (|T3|) has been obtained using dimensional regularisation and minimal 
subtraction of the ultraviolet divergence. As it is well known, only the part 

2 

proportional to log{^^) does not depend on the renormalisation scheme. 

Let us now consider the Q^i^z) term in OPE ([T^) . In general Q^i^z) will 
contain terms proportional to the VEV of the dimension-4 operator (^G^) 
and terms involving dimension-4 correlators which may arise from bifocal 
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power corrections. BPC's are taken into account to avoid ill-defined infrared 
singularities. Following this can be done by introducing the "effective 
propagator" in the form: 



'G'){k.x) A., [log 



3 ■ 36 ' TT 



+ 27i?] 



ab 



^—{k.x)^,Aflmi- 



4807r3 



+ 



(14) 



where eclipses stand for terms involving higher powers of {k ■ x) and k"^, and 
m|. = 1 GeV^, sji= 1.4 GeV^ are parameters found in This leads to the 
coefficient function Qj{z) 



Ana ,a. 



, i-G') E 4 ^ /' cos(« ■ z) [56iu) - 6iu) + 16[-]+ - 4 
144 TV ~r ^ Jo L 



u 



q y 71 



sr 

_p2 



2lE 



+ 



1207r3 



(15) 



where for any function f(u) 



du 







' du^{f{u)-f{Q)) 



(16) 



In the case of a logitudinally polarized target photon the calculation is much 
less complicated because infrared divergences are absent. The dimension zero 
coefficient Qoi^) is given by 



12a 



du sm{u ■ z) uu , 



(17) 



while for the dimension four coefficient Q^{z) we obtain 



'??! 



cos t! ■ Z) 



(18) 



At this point it is possible to compare our exact results with that of [@], 
where the photon structure function F2{u) was only calculated in the region 



6 



of intermediate u. Within the logic of the present calculation F2{u) can be 
obtained from the expression for loffe-time distribution by extracting the 
integrand in and expanding it to the first order in z. As can be easily 
seen our dimension zero coefficients Qq{z) and Qq{z) are consistent with the 
results of ref. 0. For comparison of the dimension four coefficients we have 
to neglect all singular contributions to F2{u) concentrated at the boundary 
values u = and u = 1. In general such singular contributions are crucial in 



calculations of moments of the structure function 0, |T0|. Certainly, in the 
case of loffe-time distributions they cannot be neglected. In ref. 0] BPC's 
were not taken into account. However, the dimension four contribution can 

be reproduced if the [^]+ distribution in (|15D is interpreted simply as function 

i " 

u ' 

To continue our predictions for transversally polarised on-shell photons 
we construct a Mandelstam dispersion representation in (— p^): 



oo 



Q'u) = M.) + / + r r . -"''y"'',. . (19) 



/o S ~ Jo Jo {Si — p'^){s2 — P^) 

In the region sr < —p^ <^ y?, spectral densities Ai{z;s), A2{si,S2,z) can 
be found through discontinuities, in p\ = p"^ and = + of loop 
diagrams which contribute to M^^y{p, k), equation (||). We make the standard 
approximation to represent ^40, Ai and A2 by a contribution from the p meson 
plus a continuum contribution with a treshold Sq. In this way we obtain 

Aoiz) = Coiz) + C,iz) 

Aiis;z) = 

A2{si,S2;z) = 6{si - ml)6{s2 - ml) (^C2iz) + -slCi{z)J 

+ S{si- S2)Q{tJ-'^ - si)e{si- So)y/s^Ci{z) 

- Cs{z)6{si-ml)Qi^i^-S2)e{s2-SR)-, (20) 



S2 



with the coefficient functions 



Cq{z) = — — V'e^ / du siB.{u ■ z) [{u'^ + u^) log [u u) + 1] 
q JO 

Ci(z) = — [ du sm(u ■ z) iu^ + u^) 

n ^ Jo 
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C2{Z) 



Area 
"Til 



IT 



du cos(tt ■ z) 5S{u) — S{u) + 16[ 



5. 
3' 



+ 



R 



1207r3 



e„ - z . 
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C3(^) = 47ra(^G2)^ 
This provides a representation of ""^ 



(21) 



z, —p ) valid for all —p 



C2{z) + lslCr{z) 



p 



2^2 



+ 



Co{z) + C^{z) + 
1 



log 
log( 



Si?^ - p 



+ 



p^ 



So -p^ 



Cliz) 



Csiz) 



(22) 



(-p2) (m2 - p2) ' 

which can be used to extrapolate to = 0. For numerical calculations we 
chose the normalisation scale fi"^ = 4 GeV^, and the other parameters equal 
to their standard values sq = 1.5 GeV^, sr = 1.4 GeV^, m|. = 1 GeV^, 
{^G^) = 0.012 GeV^ 

Note that the representation ( p2|) has a clear physical interpretation. The 
first two lines arise from diagonal p meson to p meson and continuum to 
continuum transitions. The contribution in the third line is due to an off- 
diagonal p — > continuum transition. The p ^ p transition is responsible 
for only half of the distribution at small z, the rest coming from continuum 
—>■ continuum and, much smaller, non-diagonal transitions. For small z, in 
the region where the present calculation can be trusted, the model loffe-time 
distribution follows practically a straight line i.e., it is practically determined 
in this region by the first derivative at 2; = which equals to the momentum 
fraction carried by quarks. 

The next step is to supplement the model loffe-time distribution for an 
on-shell photon by a pointlike contribution according to Witten's analysis, 
see []1|. Within the present formalism it amounts to the replacement of the 
factorisation scale /x^ by the virtuality of the hard photon Q^, and adjusting 
the —p"^ independent term Co{z) accordingly 



Coiz) 



3a 



du sm{u ■ z) u [{u^ + v^) 



1 



log(— ) + - 2] (23) 
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which reproduces exactly the imaginary part of the box graph 0. Further- 
more, for comparison with relatively low = 5.2 GeV^ data it is appropriate 
to neglect the perturbative evolution and treat the resulting loffe-time distri- 
bution as a low-scale, non-perturbative input. As can be seen from Figure 4 
the structure function of the real photon turns out to be larger in the region 
of small z than parametrisations |jl2| of the available F2 data. 

As a consequence the momentum fraction carried by quarks inside the 
transversal polarised p meson comes out too large as well. The loffe-time 
distribution Q'''^{z) for quarks in the transversal polarised p meson is related 
to the double p meson pole, renomalised by the coupling constant of the 
7 — > p transition 

C2{z) + \ sl C,{z) = Y.4^-^K Q'^'i^) , (24) 

^ q 9V 

where ^ ~ 1.27. The model thus predicts that essentially all longitudinal 
momentum is carried by quarks. Similiar problems are known from QCD 
sum rules calculations of the magnitude of longitudinal momentum fraction 
carried by quarks in a hadronic target 0, , although the recent calculation 
of the gluonic momentum fraction gave a reasonable result []TT|. Moreover, 



due to universality of OPE, in each case the problem can be traced down 
to the apparently too small magnitude of negative bilocal correction (|T^). 
Clearly, we have ecountered the same situation in the present context, where 
it can be equivalently related to a too small negative contribution from the 
non-diagonal transitions. When we ad hoc renormalize our loffe-time distri- 
bution by a factor 0.6 - 0.7, the plausible value of the momentum fraction 
carried by quarks in hadronic targets, the result agrees nicely with the F2 
data, leaving still some room for sea quark contributions. Of course such a 
step cannot be justified within the present formalism. 
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Figure captions 

Fig. 1 Real photon structure function -F2{z) (tliick line) as a function of 



the loffe-time z at = 5.2 GeV^ in comparison with a fit |jT2| of 
the available F2 data (dashed line). The dotted lines indicate the 
uncertainty of the F2 data. 
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Figure 1 
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